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ABSTRACT 

To  develop  an  analytical  method  for  quantifying  the  growth  behaviour  of  short  cracks 
embedded  in  notch  plastic  zones  for  power  law  strain  hardening  materials,  a  strain-based 
cohesive  zone  model  is  proposed  in  which  the  conventional  equilibrium  equation  in  the 
stress-based  model  is  replaced  by  strain  compatibility.  A  Comparison  with  finite  element 
results  shows  that  this  strain-based  model  provides  accurate  values  of  the  crack-tip-opening 
displacement  for  applied  strains  up  to  four  times  the  yield  strain  under  general  yielding 
conditions.  Furthermore,  it  is  shown  that  the  cohesive  stress  determined  by  a  method 
proposed  in  this  work  gives  better  results  than  the  existing  method,  which  are  appropriate 
only  for  small-scale  yielding  conditions. 
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A  Strain-based  Cohesive  Zone  Model  for  a 
Crack  in  a  Power-Law  Material  under  Grossly 

Plastic  Conditions 


Executive  Summary 

Modelling  of  the  growth  behaviour  of  a  short  fatigue  crack  at  a  notch  root,  with  the 
crack  being  engulfed  by  the  notch  plastic  zone,  remains  a  challenging  issue,  especially 
for  power-law  materials.  Several  schematic  models  have  been  put  forward  to  explain 
and  quantify  the  notch  fatigue  crack  growth  behaviour,  including  the  effect  of  notch 
plasticity  and  the  lack  of  plasticity-induced  crack  closure.  Attempts  have  been  made  to 
explore  the  applicability  of  elastic-plastic  fracture  mechanics  parameters,  including  the 
/-integral,  strain-intensity  factor,  etc.,  showing  only  limited  success. 

For  the  case  of  elastic-perfectly  plastic  materials,  a  strain-based  implementation  of  the 
Dugdale  model  has  recently  been  developed  for  applications  involving  a  small  crack 
engulfed  within  a  notch  plastic  zone,  but  this  work  is  limited  to  the  elastic-perfectly 
plastic  material.  In  the  present  work,  the  strain-based  Dugdale  model  is  further 
extended  to  account  for  the  effect  of  materials  strain  hardening.  A  comparison  with 
finite  element  results  shows  that  this  strain-based  model  provides  accurate  values  for 
the  crack-tip-opening  displacement  for  applied  strains  up  to  four  times  the  yield  strain 
under  generally  yielding  conditions.  Furthermore,  a  new  method  is  proposed  to 
determine  the  cohesive  stress  pertaining  to  strain  hardening  materials,  which  has  been 
shown  to  provide  better  results  than  the  existing  methods,  which  are  appropriate  only 
for  small-scale  yielding  condition. 

Results  of  the  present  research  provide  a  practical  and  physically  plausible  approach 
for  extending  the  scope  of  current  predictive  software  for  fatigue  crack  growth  based 
on  the  Dugdale  model  to  power-law  strain  hardening  materials,  including  conditions 
of  large-scale  yielding.  This  new  method  will  enable  a  more  accurate  and  reliable 
prediction  of  crack  growth  life,  which  provides  the  basis  for  setting  the  inspection 
intervals  of  aircraft. 
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NOMENCLATURE 


a  =  crack  length 

b  =  any  point  between  the  original  crack  tip  a  and  current  crack  tip  c 
c  =  length  of  fictitious  crack  c =a+r,, 
dn  =  displacement  variable 
E  =  Young's  modulus 

E  =  generalised  Young's  modulus  =Efl- tr  for  plain  strain,  E  for  plain  stress 
Ex  =  secant  modulus 

G(x,c )  =  Green  function 
h  =  half  height  of  the  plate 
k  =  strain  gradient 
Ka  ,  Ka<  =  stress-intensity  factors 
L  =  elastic-plastic  boundary 
rp  =  plastic  zone  size 

(r,6)=  polar  coordinates  centred  at  the  crack  tip 

uv  =  displacement  in  y  direction 

w  =  width  of  the  plate 

a  =  the  plastic  constraint  factor 

8m  =  crack  tip  opening  displacement 

emax  =  maximum  applied  strain 

f ,  (J'j  =  strain  and  stress 

£ep'Geq  ~  equivalent  stress  and  strain 

c 7y  =  uniaxial  yield  stress 

CT  ,  cr0  =  hypothetical  stress  distributions 

<7n  =  necking  stress 

<T0  =  cohesive  stress 

v  =  Poisson  ratio 

Vep  =  the  effective  Poisson's  ratio 
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1.  Introduction 


Modelling  of  the  growth  behaviour  of  a  short  fatigue  crack  at  a  notch  root,  with  the 
crack  being  engulfed  by  the  notch  plastic  zone,  remains  a  challenging  issue  (Suresh, 
1991;  Shin,  1994).  There  are  two  major  technical  issues:  (i)  to  develop  a  suitable 
correlating  parameter  for  the  crack  growth  rate,  especially  when  the  entire  crack  is 
embedded  within  grossly  yielded  material  (Wang  and  Rose,  1999),  and  (ii)  to 
develop  computationally  efficient  and  robust  methods  for  computing  this  correlating 
parameter,  taking  into  account  notch  plasticity.  Several  schematic  models  have  been 
put  forward  to  explain  and  quantify  the  notch  fatigue  crack  growth  behaviour, 
including  a  postulate  the  effect  of  notch  plasticity  (Smith  and  Miller,  1978)  and  the 
lack  of  plasticity-induced  crack  closure  (Newman,  1982).  Attempts  have  been  made 
to  explore  the  applicability  of  elastic-plastic  fracture  mechanics  parameters, 
including  the  /-integral,  strain-intensity  factor,  etc.,  showing  only  limited  success 
(Shin,  1994).  For  stage  II  crack  growth  in  a  ductile  material,  the  cyclic  crack-tip 
opening  displacement  (ACTOD)  has  been  shown  to  be  a  promising  correlating 
parameter  (Wang  and  Rose,  1999),  capable  of  bridging  the  gap  between  mechanically 
small  cracks  (Wang  and  Rose,  1999)  and  long  cracks  (Guo  et  al,  1999).  Here  the  usual 
operational  definition  for  CTOD  (Shih,  1981)  is  adopted.  However,  computation  of 
the  crack-tip  opening  displacement  for  a  crack  embedded  in  a  notch  plastic  zone 
represents  a  significant  challenge  as  detailed  finite  element  calculations  using 
extremely  fine  mesh  and  incremental  plasticity  theories  are  required.  Therefore  an 
analytical  method  is  required  to  calculate  the  crack  tip  opening  displacement. 

Although  the  finite  element  method  could  be  used  for  elastic-plastic  stress  analysis  of 
growing  fatigue  cracks,  thereby  providing  a  basis  for  correlating  and  predicting 
crack  growth  rates  for  arbitrary  geometry  and  loading  conditions,  the  routine 
application  of  this  approach  is  too  computationally  intensive  and  fraught  with 
difficulties  in  detailed  implementation  to  be  a  practical  option.  Instead,  a  more 
promising  alternative  is  to  use  FE  results  to  evaluate  and  to  guide  the  interpretation 
or  implementation  of  analytical  models  such  as  the  Dugdale  model  and  its  variations 
(Dugdale,  1960;  Rich  and  Roberts,  1968;  Gowda,  1972;  and  Haddad,  et  al,  1978)  and 
thereby  to  achieve  a  more  computationally  efficient  predictive  tool.  Indeed, 
considerable  success  has  been  achieved  already  with  this  approach,  most  notably 
through  the  work  of  Newman  (1982,  1992,  1998)  and  Galatolo  (1996).  However,  it  is 
recognised  by  Newman  (1998)  that  his  current  implementation  of  the  Dugdale  model 
in  FASTRAN  II  does  not  cater  for  cases  of  extensive  notch  plasticity.  This  limitation 
provided  a  major  motivation  for  the  present  work. 

In  previous  work  (Wang,  Chen  and  Rose  1999,  2002a),  the  stress-based  Dugdale 
models  were  shown  to  predict  values  for  the  plastic  zone  size  rp  and  the  crack  tip 

opening  displacement  8M  (CTOD)  which  agree  closely  with  FE  results  for  low  to 
moderate  strain  levels.  However,  predictions  diverged  from  the  FE  results  at  higher 
strain  levels  when  gross  yielding  occurs.  Furthermore,  the  strain  intensity  factor 
method  was  found  to  give  slightly  better  estimates  of  CTOD  than  the  stress-based 
Dugdale  models,  but  significantly  underestimate  the  plastic  zone  size. 


A  strain-based  Dugdale  model  (Wang,  Rose  and  Chen  2002b)  was  shown  to  give 
accurate  values  for  both  rp  and  8M  for  strain  levels  up  to  four  times  yield  strain. 
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with  no  evidence  of  significantly  diminishing  accuracy.  This  model  can  be  applied 
equally  for  plane  strain  as  for  plane  stress  with  a  constraint  factor.  However  the 
analysis  was  only  limited  to  elastic-perfectly  plastic  material  behaviour. 

The  aim  of  this  report  is  to  present  a  rational  basis  for  extending  this  strain-based 
Dugdale  model  to  the  power-law  material  response  under  a  specified  remote  strain. 
In  order  to  simulate  the  crack  behaviour  at  a  notch  root,  the  influence  of  a  strain 
gradient  become  more  important,  in  addition  to  strain-controlled  remote  loading. 
Accordingly,  the  simple  case  of  an  edge  crack  with  different  edge  crack  lengths 
subjected  to  a  remote  strain  varying  linearly  along  the  crack  path  is  considered  here. 
The  accuracy  and  range  of  applicability  of  this  strain-based  Dugdale  model  is 
assessed  by  direct  comparison  with  finite-element  (FE)  computations  detailed  in 
Section  2.  An  analysis  of  the  prospective  stress  in  the  uncracked  plate  is  presented  in 
Section  3,  while  Section  4  presents  an  extension  of  the  Strain-Based  Dugdale  model. 

2.  Finite  Element  Analysis 


There  are  two  important  features  which  need  to  be  accounted  for  in  modelling  the 
plastic  deformation  associated  with  a  small  crack  at  a  notch  root:  (i)  a  notch  strain 
field  which  varies  with  distance  from  the  notch  root,  and  (ii)  the  effect  of  the  free 
surface  at  the  notch  root  on  the  in-plane  constraint  prevailing  within  the  crack  tip 
plastic  zone. 


To  focus  on  these  two  features  in  their  simplest  form,  consider  an  edge-cracked  plate 
subjected  to  a  linearly  varying  remote  strain.  The  cracked  plate  has  the  following 
dimension:  hi  w  =  1 ,  and  at  w  —  0.05 ,  as  shown  in  Fig.  1,  where  the  crack  length  is 
denoted  as  a,  the  width  of  the  plate  and  the  half  of  height  of  the  plate  as  w  and  h, 
respectively.  To  simulate  a  linearly  varying  strain  distribution,  the  top  and  bottom 
edges  of  the  plate  are  displaced  according  to, 

uy{x)  =  hem3X{\  +  kx)  (1) 

where  ^denotes  the  maximum  strain  and  k  denotes  the  normalised  strain 
gradient: 

k  =  7td£yyW/dx  (2) 

The  remote  applied  strain  is  thus, 

M*)  =£max  (!  +  **)  (3) 

In  this  work,  the  minimum  strain  at  the  right  side  is  one  fourth  that  of  maximum 
strain  and  opposite  in  sign.  Therefore,  the  strain  gradient  is  k  —  -5/4 w . 


The  computation  was  carried  out  using  a  finite-element  code  (ABAQUS,  1998).  The 
material  is  assumed  to  obey  a  power-law  strain  hardening: 


G  eg  y 

cr. 


(4) 


with  Poisson  ratio  v  =  0.3,  and  (JY  IE  =  0.0005 .  Here  Gr  denotes  the  uniaxial  yield 
stress,  Geq  £eg  denote  equivalent  stress  and  strain,  and  a  and  n  denote  the  material 


parameters.  In  the  present  work,  various  strain  hardening  components  will  be 
considered  (n= 3,  6,  and  10)  with  a  fixed  hardening  coefficient  a  =  0.01.  The 
deformation  was  assumed  to  be  two-dimensional  plane  strain.  Due  to  symmetry. 
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only  half  plate  was  analysed.  Eight-noded  quadrilateral  elements  were  employed, 
and  the  mesh  is  shown  in  Fig.  2.  The  mesh  near  the  crack  tip  was  refined  to  capture 
the  crack-tip  plastic  blunting  behaviour. 

The  problem  depicted  in  Fig.  1  has  been  analysed  for  power-law  hardening  materials 
(Al-Ani  and  Hancock,  1991),  where  a  unique  asymptotic  field  exists  at  the  crack  tip.  It 
was  found  the  crack-tip  deformation  is  no  longer  dominated  by  the  HRR  field  when 
plastic  yielding  engulfs  the  entire  crack.  This  implies  that  the  conventional  elastic- 
plastic  fracture  mechanics  approaches  based  on  /-integral  may  fail  to  provide  an 
adequate  account  of  the  near-tip  plastic  deformation.  Therefore  a  new  method  is 
called  for  to  provide  an  accurate  assessment  of  the  plastic  deformation  at  the  tip  of 
small  crack  embedded  in  a  fully  plastic  strain  field. 

For  the  shallow  crack  shown  in  Fig.  2,  as  the  applied  maximum  strain  increases, 
plastic  deformation  starts  to  develop  in  two  separate  regions,  ahead  of  the  crack  tip 
and  away  from  the  crack,  see  Fig  3(a),  (b).  Finally,  when  a  certain  strain  is  reached, 
the  two  plastic  regions  merge  into  one,  as  shown  in  Fig.  3(c). 

A  measure  of  the  crack-tip  plastic  deformation  that  has  important  implications  for 
fatigue  crack  growth  is  the  crack-tip  opening  displacement.  The  parameters  of 
primary  interest  in  the  present  work,  for  the  purposes  of  comparison  with  the 
analytical  estimates  based  on  the  Dugdale  model,  are  the  plastic  zone  size  rp  and  the 

crack-tip  opening  displacement  dM  .  Numerical  values  for  these  parameters  were 
derived  from  the  FE  results  by  the  procedures  shown  below.  The  crack-tip  opening 
displacement  is  defined  as  the  crack  face  separation  where  the  45-degree  lines  drawn 
back  from  the  crack-tip  intersect  the  deformed  crack  face  (Shih,  1981),  as  illustrated  in 
Fig.  4,  which  can  be  mathematically  expressed  as 

SM=2uy(r*,0=Ji)  (5) 

where  r  *  is  the  root  of  the  following  equation, 

r-ux(r,d  =7t)  =  uy(r,0  =k)  (6) 

where  (r,6)  denote  polar  coordinates  centred  at  the  crack  tip.  As  shown  by  Shih 
(1981),  this  construction  gives  a  CTOD  that  matches  the  Dugdale  model  for  the  small- 
scale  yielding  case  under  plane  stress  conditions. 
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(a)  (b) 


Figure  1  An  edge-cracked  plate  subjected  to  a  linearly  varying  strain  field  showing  (a)  the 
strain  distribution  and  (b)  the  attendant  stress  distribution  for  power-law  strain 
hardening  material. 


Figure  2  Finite  element  model  of  half  plate 
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Figure  4  The  definition  of  crack-tip  opening  displacement 


3.  Stress  distribution  in  the  un-cracked  plate 


For  the  nominal  elastic  stress  distribution  in  an  uncracked  plate  subjected  to  the 
linearly  varying  strain  given  by  Eq  (3),  under  the  plane  strain  conditions,  two  issues 
need  to  be  addressed:  (i)  accounting  for  the  effects  of  stress  triaxiality;  (ii)  selecting  an 
appropriate  constraint  factor  a  for  the  limit  stress  <J0  in  the  Dugdale  model.  The 
elastic  plastic  stress  distribution  in  the  absence  of  a  crack  can  be  calculated  using 
Hencky's  deformation  theory  of  plasticity  (Chakrabarty,  1987),  this  is  because  the 
deformation  state  in  the  uncracked  plate  is  a  proportional  one:  ratios  of  stress 
components  and  ratios  of  strain  components  are  independent  of  applied  strain. 


For  the  plane  strain  case,  £_  =0,  the  non-zero  strain  components  are:  e  and  £ vv  , 
and  the  associated  non-zero  stress  components  are  <7  vr  and  <7_.  .  The  condition  of 
symmetry  implies  that  <7 xy  v  =0  from  the  equilibrium  equations  CTu  t  =-<7u  v  =  0, 
so  that  in  view  of  the  symmetry  condition  a xv{y  =  0)  =  0  and  of  the  free-surface 
boundary  condition  <7„(.r  =  0)  =  0,  one  finds  <7XX  =  Txv  =0  throughout  the  body. 
According  to  deformation  theory  of  plasticity,  the  two  non-zero  stresses  are  related, 

<*=  =  %<*,,.  (7) 

where  vep  denotes  the  effective  Poisson's  ratio, 

vep  =  0.5- (0.5 -V)^-  (8) 

Here  Es  is  the  secant  modulus,  defined  as: 


The  equivalent  stress  Geq  is,  by  definition: 


(9) 
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V,q=^-Vc„+Vl>°,< 

and  according  to  the  generalised  Hooke's  law: 

'~v„, 

Substituting  Eqs  (9, 10)  into  Eq  (11),  we  can  get: 


(10) 

(11) 

(12) 


Here,  v  needs  to  be  determined  iteratively.  Inserting  Eq  (8-12)  into  Eq  (4): 


£  it  (x)  is  the  applied  far  field  strain  from  Eq  (3)  and  v  =0.3,  effective  poison  ratio 
V  ^  can  be  readily  determined  by  solving  the  above  non-linear  equation.  Then 
g  ,<Tn.  and  (7cqcan  be  determined  from  Eq  (12),  Eq  (11)  and  Eq  (10)  respectively. 
Therefore,  the  distribution  of  the  elasto-plastic  stress  can  thus  be  expressed  as 


<MX) 


E  eyy(x)  , 


x  <  L 
x  >  L 


(14) 


where  E  denotes  the  generalised  Young's  modulus  ( =  E  for  generalised  plane 
stress,  or  £7(1 -V2)  for  plane  strain,  v  denoting  Poisson's  ratio),  L  represents  the 
elastic-plastic  boundary: 


(15) 


The  comparison  between  the  analytical  and  FE  results  is  shown  in  Fig.  5.  The  elasto- 
plastic  stress  distribution,  which  would  prevail  in  the  absence  of  a  crack,  is 
independent  of  y  for  plane  strain  conditions. 
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Figure  5  A  comparison  of  prospective  stress  between  the  predictions  and  FE  results.  Lines 
denote  predictions,  and  symbols  FE  results. 


4.  Strain-based  Dugdale  Models 

The  stress-based  Dugdale  model  has  been  evaluated  in  (Wang,  Chen  and  Rose, 
2002a),  and  it  was  found  that  the  standard  implementation  of  this  model  is  unable  to 
adequately  characterise  the  crack-tip  opening  displacement  dM  and  the  plastic  zone 
size  rp.  Consequently,  a  strain-based  Dugdale  approach  was  developed  (Wang, 
Rose  and  Chen,  2002b)  for  a  crack  in  a  perfectly  plastic  material. 

4.1  Stress-Based  Dugdale  Model 

For  a  plate  subjected  to  a  combined  bending  and  tension,  with  the  prospective  stress 
having  been  calculated,  the  crack-tip  opening  displacement  can  now  be  calculated 
using  the  following  method. 

According  to  Dugdale's  cohesive  zone  concept  the  plastic  zone  size  rP  of  a  crack 
length  a  is  given  by  the  root  of  the  following  equation. 
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*„</>>+ *<,„<'V>  =  0  <16> 

If  a  represents  the  original  crack  size,  and  the  rf,  represents  the  plastic  zone  size  in 
front  of  the  crack  tip,  and  then  c-a  +  rp  .  Eq  (16)  can  be  readily  solved  by  iteration, 

given  the  stress  distribution  along  the  crack  o(x),  here  Ka  and  Ka  are  the  stress- 
intensity  factors  for  the  two  sub-problems,  which  can  be  written  as: 


and 


Ka  (rp )  =  j‘o(x)G(x,c)dx 

1 

Ka(rp)  =  -a\'(J0G(x,c)dx 


(17) 


(18) 


where,  a  is  the  plastic  constraint  factor:  a  =  1  for  plane  stress  and  O'  =  1 .73  (Erwin, 

1960)  for  plane  strain.  The  function  G  denotes  the  stress-intensity  factor  of  a  crack  of 

half-length  c  subjected  to  a  pair  of  point  forces  acting  at  x.  The  Green  function  for  an 

edge  crack  in  a  half-plane  is  (Tada,  1985), 

~  2  1.3 -0.3  (x/cf4 

G(x,c)  = 


\2  il  /2 


(19) 


I7TC  [1-(X/C)2]' 

The  total  CTOD  is  the  sum  of  the  displacement  Ay  due  to  the  traction  <7n.  (x)  and  the 
displacement  ACT  due  to  the  cohesive  stress  cr0 : 


A  =  A„+A„,  (W) 

Where  (Tada  et  al,  1985): 

A  =  —  f  f  a  yv  (x)G(x,  a)G(a,  b)dxdb  (20) 

°  p  Ja  JO  - - 

where  b  represents  any  point  between  the  original  crack  tip  a  and  current  crack  tip  c, 
and 

=  -a  f  f  <J0G(x,b)G(a,b)cbcdb  (22) 

0  Ja  Ja 


4.2  Strain-Based  Dugdale  Model 

Like  the  Dugdale  model,  the  BCS  model  (Bilby  et  al,  1963;  Bilby  and  Eshelby,  1968) 
also  uses  a  continuous  distribution  of  dislocations  to  model  both  a  crack  and  its 
associated  plastic  zone.  Although  the  same  mathematical  formalism  applies  to  both 
models,  the  physical  interpretation  of  the  BCS  model  suggests  a  plausible  extension 
of  the  model  for  cracks  under  fully  plastic  loading.  The  dislocation  density  D(x)  in 
the  BCS  model  is  determined  from  the  following  singular  integral  equation. 


(23) 


fc  [-<7  (x)/ E‘  (0<x<tf) 

f  D(t)Q(x,t)dt  =  \  -  .  , 

J  o  [[— cr  vv  (x)  +  <70  )]/E  (a  <  x  <  c) 

where  the  kernel  Q(x,t)  has  a  Cauchy-type  singularity  for  x  =  t ,  the  model  is 

restricted  to  cases  where  the  applied  stress  <J)y  (x)  is  within  the  elastic  range.  If  we 

view  the  right  hand  side  of  Eq  (23)  in  terms  of  the  applied  strain  <7 yv  (x)  /  E  ,  rather 
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than  the  applied  stress,  which  controls  the  dislocation  density,  the  model  could  still 
be  applied  to  cases  where  the  applied  stress  exceeds  the  yield  stress,  which  would 
prevail  in  the  absence  of  a  crack.  Similarly,  the  Dugdale  model  can  be  rewritten  in 
terms  of  total  strain  in  the  same  form  as  Eq  (23)  by  defining  hypothetical  stress 
distributions  <7  and  a'0  as  follows. 


jjDU)Qx,t)dt= 

where 


J-<7  (x)/  E 
{[-<J(x)+ct0(.y)]/E 


(0  <  x  <  ct) 
( a<x<c ) 


(24a) 


a  (x)  =  E  £,/;(x)  (24b) 

<70(x)  =  c j0E  /EJx)  (24c) 

Equations  (24a-c)  are  almost  identical  to  those  for  the  stress-based  model  provided 
that  <7  (x)  given  by  Eq  (24b)  is  regarded  as  defining  the  applied  stress  which  would 
prevail  in  the  absence  of  a  crack,  and  <7'(x)  is  interpreted  as  the  new  limit  stress, 
which  is  now  a  function  of  position,  through  Eq  (24c).  Nevertheless,  the  extended- 
crack  length  c  can  still  be  determined  iteratively  from  the  condition, 

K(c)  =  0  (25a) 

with  <7n.(x, y  =  0)  now  taken  to  be  specified  by 


C7  n  (x,  v  =  0)  = 


j-cr*(jc) 

[— ct^Hct^x) 


(0  <  x  <  a) 
(a  <  x  <  c ) 


(25b) 


The  advantage  of  this  formulation  is  that  convenient  analytical  approximations  for 
the  weight  function  G(x,c )  in  Eq  24(a)  are  now  available  for  a  wide  range  of  crack 
geometrises  (Tada,  1985),  even  though  exact  analytical  solutions  are  available  for 
even  fewer  geometrise  than  for  the  dislocation  kernel.  For  example,  an  exact 
analytical  solution  for  G(x,c)  is  not  available  for  the  case  of  a  crack  emanating  from 
a  circular  hole,  or  even  for  the  limiting  case  of  an  edge  crack  in  a  half-plane. 
However,  for  the  latter  case,  Tada  (1985)  gives  the  convenient  approximation  of 
acceptable  accuracy  shown  in  Eq  (19). 


Given  an  explicit  representation  such  as  Eq  (19)  for  G(x,c) ,  one  can  derive  the 
extended-crack  length  c ,  and  hence  the  plastic-zone  size,  r  =c  —  a,  for  the  strain- 

based  model.  The  crack  opening  displacement  <5(x)is  related  to  the  dislocation 
density  through, 

(*)  =  ",  (*.  y  ->  0+)  -  uy  (x,  y  ->  0-),  (26a) 

=  \iD(t)dt  (26b) 

which  can  be  determined  using  the  method  outlined  in  Section  4.1 . 


4.3  Cohesive  Stress 

The  application  of  the  Ehigdale  model  for  plane  strain  requires  the  use  of  semi- 
empirical  constraint  factors  that  have  been  extensively  investigated  and  discussed  in 
previous  work  (Newman,  1982;  1998;  Guo  et  al,  1999).  A  further  discussion  regarding 
the  appropriate  choice  of  constraint  factor  can  be  found  in  the  work  (Wang,  Rose  and 
Chen,  2002b). 
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The  last  remaining  problem  for  the  above  strain-based  Dugdale  model  of  power  law 
material  is  to  determine  the  cohesive  stress.  Denoting  cr0  as  the  cohesive  stress,  one 
can  simplify  the  stress-strain  curve  of  the  power  law  materials  as  elastic-perfectly 
plastic  in  which  the  yielding  stress  is  assumed  as  the  cohesive  stress.  This  cohesive 
stress  can  be  selected  by  using  three  methods:  (1)  Necking  Stress  Method:  the 
cohesive  stress  equals  to  the  necking  stress;  (2)  Equal-Area  Method:  the  cohesive 
stress  is  selected  to  give  the  same  area  below  the  stress-strain  curve;  and  (3)  HRR 
Method:  the  cohesive  stress  is  selected  according  to  a  small-scale  yielding  (SSY) 
solution.  These  three  methods  are  summarised  in  the  following  section. 


Denoting  eyas  necking  stress  (Chakrabarty,  1987),  which  is  the  root  of  equation 

—  =  0 ,  the  corresponding  strain  at  which  is  £  =  1  /  n  .  Substituting  this  solution  into 
d£ 

the  Eq  (4),  one  can  get: 

_L  =  £jl  +  J^l>)  (27) 

£yn  Gy  [Gy  ^ 

Solving  the  above  equation,  we  can  obtain  the  necking  stress  Gn  for  a  given  strain. 
Simply  we  can  directly  use  the  necking  stress  as  the  cohesive  stress,  which  is  referred 
to  as  Necking  Stress  Method  here.  Rice  [1968]  proposed  to  use  Gn  in  conjunction  with 


Dugdale  model  solution: 


(28) 


If  a  solid  curve  denotes  the  original  stress-strain  curve  of  power  law  material  shown 
in  Fig  (6),  and  a  doted  curve  denotes  a  new  stress-strain  curve  derived  from  the 
original  stress-strain  curve  according  to  the  Equal-Area  Method,  which  means  that 
the  area  below  the  solid  curve  should  be  equal  to  the  area  below  the  dotted  line.  As 
can  be  seen  in  Fig.  (6),  the  new  dotted  curve  represents  an  elastic-perfectly  plastic 
curve  in  which  the  yielding  stress  is  cohesive  stress  G 0 .  Because  the  area  under  the 
solid  curve  equals  to  the  area  under  the  dotted  curve,  we  get: 

^ Eg, 2  +  G0  (£max  — £0)  =  ~Eoy2  +  £“  <7  d£  (29) 

Here  the  maximum  stress  is  assumed  as  the  necking  stress.  Substituting  the 
corresponding  necking  strain  into  the  above  equation: 

^£((702  — O  +  C70(--£0)=  f  n°d£  (30> 

2  n  Je> 

From  the  above  equation,  the  cohesive  stress  can  be  solved.  This  method  is  suitable 
for  large  scale  yielding  (LSY)  condition,  and  we  refer  to  it  as  the  equal-area  method. 


Another  method  of  obtaining  the  cohesive  stress  is  to  use 
crack-tip  opening  displacement  CTOD  is  given  by: 


the  SSY  solution.  The 
(31) 


and  the  relationship  between  J  and  K  is  given  by 
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therefore,  the  final  formulation  of  CTOD  can  be  written  by: 
,  K2 


8\i 


where  <7 ,  is  the  yielding  stress,  and  d„  is  the  displacement  variable  which  is  a 

function  of  n  and  ) .  Recalling  the  elastic  theorem,  the  formulation  of  CTOD  is- 
E 

<34> 

Matching  the  form  of  Eq  (34)  of  the  elastic  theorem,  the  plastic  Eq  (33)  can  also  be 
rewritten  as  : 

c  K 2 

(35) 


Comparing  Eq  (33)  and  Eq  (35),  we  can  get  the  following  formulation  of  cohesive 
stress: 


(36) 

The  values  of  d„  can  be  found  in  Anderson  (1985).  This  method  is  suitable  for  small 
scale  yielding  cases  (SSY),  and  we  call  it  HRR  method.  There  are  some  differences 
between  the  results  of  LSY  and  SSY  methods  as  shown  in  Fig.  (7). 
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Figure  7  Acomparison  of  cohesive  stress  by  using  the  HRR  Method,  Necking  Stress 
Method  and  Equal- Area  Method. 


0  12  3  4 

Maximum  applied  strain  e^/Ey 

Figure  8  Comparison  betzveen  predictions  of  strain-based  Dugdale  model  and  finite  element 
results  for  crack-tip  opening  displacement  under  plane  strain  conditions.  The  solid 
lines  indicate  predictions  using  the  cohesive  stresses  obtained  from  the  Equal-Area 
method  and  the  dashed  lines  indicate  predictions  using  cohesive  stresses  from  the 
Necking  Stress  Method, and  symbols  indicate  FE  results. 
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5.  Results  and  Discussion 

In  this  section,  the  predictive  accuracy  of  the  analytical  models  presented  in  Section  2 
will  be  assessed  against  the  computational  (FE)  results.  Some  issues  will  also  be 
briefly  discussed,  along  with  recommendations  on  promising  directions  for  future 
work. 

The  analytical  estimates  of  crack-tip-opening  displacement  based  on  strain-based 
model  compared  with  FE  results  are  discussed  here.  The  controlling  parameter  is  the 
maximum  applied  strain  emax ,  normalised  by  the  uniaxial  yield  strain  £,  =cr)  IE. 
In  order  to  use  the  two-dimensional  Dugdale  model  to  characterise  the  complex 
three-dimensional  plastic  deformation  of  plane  strain  cracks,  the  limiting  stress  is 
chosen  to  represent  a  thickness-average  of  the  normal  stresses  in  the  plastic  zone. 
However,  it  has  been  reported  that  the  Dugdale  model  with  the  stress-based 
constraint  factor  provides  a  reasonable  estimate  of  the  plastic  zone  size,  but  over¬ 
estimates  the  crack-tip  opening  displacement  (Newman  et  al,  1995).  Under  large- 
scale  yielding  and  gross-section  yielding  conditions,  the  significant  under-estimation 
of  the  plastic  zone  size  by  the  stress-based  method  indicates  that,  due  to  the  loss  of 
plastic  constraint,  a  lower  constraint  factor  would  seem  more  appropriate  (Wang 
2000).  Fig.  8  shows  the  results  obtained  from  the  strain-based  Dugdale  model  of 
Section  3.  The  solid  lines  indicate  predictions  using  the  cohesive  stresses  obtained 
from  the  Equal-Area  Method  and  the  dashed  lines  indicate  the  predictions  of  the 
Necking  Stress  Method.  It  can  be  seen  that  former  predictions  are  in  better  agreement 
with  the  FE  results  than  the  latter  predictions.  This  confirms  that  the  Equal-Area 
method  of  cohesive  stress  proposed  in  this  work  pertain  to  the  problem  of  a  small 
crack  engulfed  within  a  notch  plastic  zone  for  power-law  strain  hardening  materials, 
rather  than  the  usual  choice  of  cohesive  stress  for  small-scale  yielding  conditions  at  a 
crack  tip. 


Since  the  objective  for  the  present  work  is  to  develop  a  practical  approach  for 
predicting  fatigue  crack  growth  rates  for  power-law  materials  at  notch  roots,  when 
the  crack  is  still  small  compared  with  the  notch  plastic  zone,  the  significance  of  the 
strain-based  Dugdale  model  in  this  context  and  the  requirements  for  implementation 
of  the  present  approach  in  practice  will  now  be  discussed,  along  with 
recommendations  for  promising  directions  for  future  work. 

Two  key  parameters  are  required  for  the  implementation  of  the  strain-based  Dugdale 
model,  viz,  (i)  the  limit  stress  <J0,  or  equivalently  the  constraint  factor  a ;  and  (ii)  the 
appropriate  secant  modulus  Es.  As  indicated  above,  for  arbitrary  notch  geometrise 
these  would  need  to  be  determined  from  an  FE  analysis,  but  it  would  also  be  able  to 
expect  that  adequate  approximations  could  be  developed  for  standard  notch 
geometrise,  as  discussed  by  Wang  et  al  (1999),  thereby  reducing  the  computational 
burden  in  routine  applications.  Another  requirement  is  the  availability  of  analytical 
expressions  or  approximations  for  the  dislocation  solution  Q(x;c)  in  Eq  (9a),  or 
equivalently  for  the  weight  function  G(x;c).  As  noted  in  Section  2,  Tada's  work 
(1985)  has  shown  that  analytical  approximations  for  G(x;c)  are  now  available  for  a 
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wide  range  of  notch  and  crack  geometrises,  so  that  practical  implementation  of  the 
present  approach  would  not  be  hindered  by  this  requirement 

This  strain-based  Dugdale  model  is  still  based  on  the  superposition  method  which 
reduces  the  elastic-plastic  problem  into  a  superposition  of  two  elastic  solutions  of  a 
fictitious  crack  subjected  to,  respectively,  a  prospective  stress  that  would  exist  in  the 
absence  of  the  crack,  and  a  cohesive  stress  equal  to  the  yield  stress  of  the  perfectly 
plastic  material.  In  this  work,  the  strain-based  Dugdale  model  is  extended  to 
consider  power-law  materials  by  simplifying  the  constitutive  equation  of  power  law 
material  into  an  effective  elastic-perfectly  plastic  constitutive  law,  offering  a  simple 
yet  efficient  mean  of  characterising  the  plastic  deformation  at  a  crack  tip.  In 
particular,  for  power-law  materials,  the  determination  of  cohesive  stress  has  a 
significant  impact  on  the  results.  As  addressed  in  Section  4,  there  are  three  methods 
to  determine  the  cohesive  stress,  (1)  the  necking  stress  method,  (2)  LSY  method  and 
(3)  SSY  method.  In  this  work  a  good  correlation  has  been  found  by  using  the 
cohesive  stress  determined  by  proposed  LSY  method,  in  which  the  maximum  stress 
is  selected  as  the  necking  stress,  thus  this  approach  is  recommended  to  the  large  scale 
yielding  cases.  The  bound  of  applicability  of  these  different  methods  is  not  very 
clear.  Therefore  further  work  in  this  direction  would  be  required. 


6.  Conclusion 

A  strain-based  implementation  of  the  Dugdale  model  has  been  extended  for  further 
applications  involving  a  small  crack  engulfed  within  a  notch  plastic  zone  of  power- 
law  strain  hardening  materials.  The  predictive  accuracy  of  this  model  has  been 
evaluated  by  comparison  with  elasto-plastic  finite  element  analysis.  To  retain  in  their 
simplest  form  of  the  key  features  of  strain  control  and  of  a  strain  gradient  that 
characterise  the  deformation  at  a  notch  root,  this  evaluation  focused  on  the  response 
of  an  edge  crack  subjected  to  a  linearly  varying  strain  distribution.  The  strain-based 
model  was  shown  to  give  accurate  values  for  dM  for  strain  levels  up  to  £max  / eY  =  4 , 
with  no  evidence  of  significantly  diminished  accuracy  with  increasing  strain, 
provided  that  the  cohesive  stress  is  determined  by  the  proposed  LSY  method  which 
is  appropriate  for  the  notch  plastic  deformation  in  the  absence  of  a  crack  for  power- 
law  hardening  materials. 
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